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Abstract. Interstellar superbubbles generated by multiple supernova explosions are common in 
star-forming galaxies. They are the most obvious manifestation of mechanical feedback, and are 
largely responsible for transferring both thermal and kinetic energy to the interstellar medium 
, from the massive star population. However, the details of this energy transfer remain surprisingly 
' murky when individual objects are studied. I will summarize what we currently know about 
candidate dominant processes on these scales. 
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1. Introduction 

The previous contributions showed a wealth of H i shells and supershells. How well- 
established is it, that these structures originate from mechanical feedback, namely, winds 
, and supernovae (SNe) from OB associations? If we are to understand how massive star 
J> ■ feedback affects galaxies on global scales, then it is essential that we understand the 
00 origin and evolution of superbubbles, which are the direct manifestation of feedback. 
• We refer not only to mechanical energy, but also to photoionization, and dispersal of 
' nucleosynthesis products from the SNe and massive stars. In this presentation, I use the 
definition introduced by You-Hua Chu in the 1980's: a "superbubble" is a shell structure 
originating from multiple stars, rather than a definition related to its size. 

The default model for superbubbles is the standard, adiabatic evolution that assumes 
' a hot (10^-1- K), shock-heated interior whose pressure drives the outer shell growth 
O I (Pikel'ner 1968; Castor et al. 1975; Dyson 1977). For this model, the shell parameters are 
^ determined by only three input parameters: the mechanical "luminosity" L, the ambient 
^ density n, and the age t, e.g., 

1: Ro,{L/nf'U''' , (1.1) 

„ ^ ^ « (i/,)2/5 ,-2/5 (1,2) 

5h ' ^® ^^^'^ must know conditions in the ambient environment, for example, the ambient 
pressure and density distribution. The above is the simplest, analytic representation of 
this model, but there is also a large body of work in hydrodynamic simulations of these 
objects over the past two decades (e.g., Mac Low et al. 1989; Bisnovatyi-Kogan et al. 
1989; Palous et al. 1990; Tenorio-Tagle et al. 1991; Tomisaka 1990; Silich et al. 1996; 
Gazol-Patino & Passot 1999; Strickland & Stevens 2000). 

What evidence do we seek, that "star formation," i.e., massive stars, are responsible for 
the ubiquitous superbubbles that we see in gas-rich galaxies? Since the adiabatic model 
is based on the shock-heated interior, we expect gas at multiple temperature phases, 
including hot. X-ray emitting gas; warm, 10* K photoionized gas; and also the neutral 
gas that we have seen. Second, we should also find that objects whose input parameters 
are known should follow the prescribed evolution above. Third, the statistical properties 
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of the superbubble populations should be consistent with the statistical properties of 
the putative parent star-forming regions. Fourth, and perhaps most obvious, we expect 
a one-to-one correspondence between the superbubbles and the massive star clusters. 

2. Multi-phase ISM 

As mentioned by Elias Brinks in his talk, we often do see the existence of gas at the 
different temperatures that are predicted by the adiabatic model. A recent example is 
the X-ray emission in N51 D (Cooper et al. 2004), along with photoionizod Ha emission. 
Superbubbles were already detected in X-rays by Einstein and ROSAT (e.g., Chu & Mac 
Low 1990; Wang & Helfand 1991). It has long been known that there are two categories 
of X-ray emission from these objects: X-ray bright, and X-ray dim (Chu et al. 1995). The 
former show X-ray emission in excess of what is expected from the adiabatic model, and 
the enhancements are likely caused by secondary SN blastwave impacts to the shell walls 
(Chu & Mac Low 1990; Oey 1996). 

H recombination emission is generally consistent with photoionization by the observed 
early-type stars, although density-bounding and shock-heating contributions are also 
often factors (Hunter et al. 1995; Oey et al. 2000). We usually sec Hq emission on the 
interior of H i shells, consistent with the model that star formation is responsible for the 
superbubbles. Kim et al. (1999) studied the Hi shells of the Large Magellanic Cloud 
(LMC) and were even able to suggest an evolutionary sequence defined by relative Ha 
and H i morphology. 

Intermediate-temperature ions are also present. C iv and Si iv, are often seen in 
absorption in the lines of sight toward massive stars within superbubbles (Chu et al. 
1994). Recently, O Vi was also detected by FUSE in the line of sight toward one LMC 
superbubble (Danforth & Blair 2006). 

3. Dynamics of Individual Superbubbles 

If the input mechanical power, ambient conditions, and age of the superbubbles can 
be determined, we can test for consistency with the adiabatic evolution. Such studies 
are possible for objects in the Milky Way and Magellanic Clouds, where the stellar 
population can be resolved, and these studies invariably show the existence of a growth- 
rate discrepancy such that the shells are much smaller than predicted by the apparent 
input parameters (Saken et al. 1992; Brown et al. 1995; Oey & Massey 1995; Hunter et 
al. 1995; Oey 1996; Cooper et al. 2004). The problem has been known for many years, 
and is even seen in single-star nebulae generated by Wolf-Rayet stars (e.g., Cappa et al. 
2001, 2005). 

What is the current status in resolving this growth-rate discrepancy? There are a num- 
ber of possible important factors. Those discussed below are the most likely candidates. 
Combinations of multiple effects may also be at play, including others, such as viscous 
drag (I. Goldman, private communication), that we do not discuss fully below. 

3.1. Input power overestimated? 
The first possibility for resolving the superbubble growth-rate discrepancy is that the 
input mechanical luminosity has been overestimated. Stellar wind mass-loss rates M are 
especially suspect, since it has long been suggested that clumping in the winds leads to 
overestimates in M from radio continuum measurements (e.g., Hillier 1991; Nugis et al. 
1998). In recent years, this appears to be confirmed by X-ray line profile fitting (Cohen 
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et al. 2006; Miller et al. 2002; Kramer et al. 2003). FuUerton et al. (2006) also find the 
same result from P V line profiles observed by FUSE. Their analysis of this dominant ion 
suggests that the overestimates in M may be as high as two orders of magnitude in some 
cases! These overestimates in M are relevant primarily to the youngest superbubbles, 
whose evolution is still dominated by stellar winds instead of SNe, which applies to most 
dynamical studies of optical or X-ray selected superbubbles. 

3.2. Ambient density underestimated? 

As seen in equations 11.11 and 11.21 an underestimate in the ambient density has an equiva- 
lent effect to an overestimate in L. Oey et al. (2002) mapped the immediate environment 
of three Ha-selected LMC superbubbles to determine whether the neutral gas environ- 
ment was unusually dense. We found an extreme range in conditions for the three objects, 
with one essentially in an H i void, another nestled amongst a number of H i clouds, and 
the third in a region with no obvious relation to the observed H i. This lack of any sys- 
tematic effect suggests that the ambient density may not be the primary source of the 
growth-rate discrepancy, although it also demonstrates that the ambient environment is 
more complex than assumed. 

We also note that if a superbubble originates in a higher-density cloud, then a mini- 
blowout from the cloud can also accelerate the shell's observed expansion velocity relative 
to its radius (Oey & Smedley 1998; Mac Low et al. 1998). Hence it is possible to reproduce 
unusual observed kinematics in particular objects. 

3.3. Ambient pressure underestimated? 

Another important environment parameter is the interstellar pressure, which counteracts 
the growth of the shells. In general, models assume a fiducial ambient P/k on the order of 
10^ cm^'' K or less. The various sources of ambient pressure can be broadly described as 
the thermal pressure, magnetic pressure, turbulent pressure, and cosmic ray pressure. We 
note that all of these correlate with star- formation rate (SFR). Likewise, superbubbles 
tend to be prominent in active star-forming galaxies like the LMC and localized active 
environments. Indeed, the very existence of star-formation in these regions implies higher 
local pressure (M. Mac Low, private communication). It is therefore plausible that the 
ambient pressure for superbubbles may have been systematically overestimated. 

Oey & Garcfa-Segura (2004) discuss this possibility and present 2-D hydrodynamic 
simulations of six LMC superbubbles for ambient P/k = 10'* and 10^ cm"'^ K. The mod- 
els are generated with the same input mechanical luminosity that is estimated from the 
observed stellar population by Oey (1996). For the lower pressure, the simulated radial 
density profiles show an extended photoionized morphology that is inconsistent with the 
observed Ha data. In contrast, the simulations assuming ambient P/k = 10^ cm^^ K 
show no extended Ha emission, and an Ha morphology that agrees well with the ob- 
servations. The observed and modeled velocity structures are similarly more consistent 
with the models for the high-pressure environment in all cases. The H i radial profile 
predictions also differ, and can be used to further test these models. 

3.4. Radiative Cooling? 

If the hot superbubble interiors are in fact losing energy by radiative cooling, then the 
growth will not keep pace with the adiabatic model. Mass-loading into the hot interior 
via evaporation from the shell walls or ablation of clouds and clumps can drive radiative 
cooling (e.g., Hartquist et al. 1986; Arthur & Henney 1996). Alternatively, an increase 
in metallicity due to the injection of nucleosynthesis products from the parent SNe can 
also enhance the cooling rate. Silich & Oey (2002) show that the X-ray luminosity of 
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a low-metallicity {Z — 0.05Zq) superbubble can be increased by almost an order of 
magnitude, simply by products from 3-4 SNe. Silich et al. (2001) examine this issue for 
starburst superwinds. However, the observed X-ray luminosities generally do not appear 
to suggest anomalous cooling from the superbubble interiors (e.g., Chu et al. 1995). 



Another factor whose importance has been underemphasized, is the transfer of superbub- 
ble energy to cosmic rays. Multi-SN superbubbles are especially efficient at accelerating 
cosmic rays because the blastwaves expand into a pre-heated environment (e.g., Parizot 
et al. 2004). The superbubble interiors thus harbor strong MHD turbulence and mag- 
netic fields (e.g., Bykov 2001; Bykov & Toptygin 1987), which are needed for cosmic ray 
acceleration. Because of the multiple SNR shocks, superbubbles also promote multiple 
accelerations, which can push cosmic rays to higher energies (Parizot et al. 2004; Klepach 
et al. 2000) . The cosmic ray energy distributions and isotope abundances are broadly con- 
sistent with superbubble origins. While the role of superbubbles in explaining cosmic ray 
properties has been recognized for some time, the effect of energy transfer on the parent 
objects themselves has only been discussed recently, and rough estimates suggest that a 
few tenths of superbubble kinetic energy could be lost to cosmic rays (Parizot et al. 2004; 
Bykov 2001). For individual young SNRs, simulations predict growth deviations > 10% 
and reduced X-ray luminosities (Ellison et al. 2004). This energy sink applies more to 
SN-dominated superbubbles, rather than stellar wind-dominated objects. 

4. Properties of global mechanical feedback 

Another approach to evaluating the adiabatic model is to compare observations and 
predictions for the statistical properties of superbubbles in galaxies, based on the known 
global star-formation properties. Oey & Clarke (1997) derived size distributions and 
expansion velocity distributions for extremes in star-formation history and star-cluster 
mass functions. The latter, which produces the H ii region luminosity function, is gener- 
ally a robust power law with a dependence of in differential form (e.g., Efremov & 
Elmegreen 1997; Oey & Clarke 1998). This yields a differential size distribution. 



We find excellent agreement with these relations for the H i shell population found by 
Staveley-Smith et al. (1997) for the SMC. More recently, Hatzidimitriou et al. (2005) 
updated the SMC shell catalog and re-examined these relations, inferring a different 
star-formation history, but one that is still consistent with a star- formation origin for the 
structures, plus adiabatic evolution. H i shell size distributions for a few other galaxies 
also have been examined, for example, the Milky Way (S. Ehlerova, these proceedings), 
LMC (Kim et al. 1999), M31, M33, Holmberg II (Oey & Clarke 1997), and NGC 2403 
(Thilker et al. 1998; Mashchenko et al. 1999). These other catalogs are generally consis- 
tent with equation 14. 11 but the statistics are more incomplete. 

Analysis of the global size distributions leads to the definition of a critical star- 
formation rate, above which the shells merge and shred the neutral ISM, generating 
pressure-driven outflows from their parent galaxy disks (Clarke & Oey 2002): 



3.5. Energy Transferred to Cosmic Rays? 



N{R) dR oc R^^ dR 
and corresponding distribution in expansion velocities, 

N{v) dv oc u""^/^ dv 



(4.1) 



(4.2) 



SFRcrit — 0.15 MisM,io cr^,io//d Mq yr ^ 



(4.3) 
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where Mism.io and av,io are the mass of the ISM in units of 10^" M© and the thermal 
velocity dispersion in units of 10 kms~^, respectively. The geometric correction factor 
fd is on the order of unity. Is this expected shredding of the neutral ISM for high SFR 
morphologically apparent in the H i datasets? 

We can compare the Hi maps and shell catalogs for the LMC (Kim et al. 1999) 
and SMC (Staveley-Smith et al. 1997). The two surveys were both carried out with the 
Australia Telescope, at similar depth and spatial resolution. Although the LMC is a 
larger galaxy and has a much higher star-formation rate than the SMC, the number 
of coherent H l shells identified in the LMC survey is only 126, roughly one-quarter of 
the 509 (Hatzidimitriou et al. 2005) found in the SMC. Ordinarily, we would expect 
a much larger number of shells in the LMC than in the SMC. Morphologically, there 
is a noticeable contrast between these two galaxies: the LMC H i has a filamentary 
appearance, consistent with shredding and compression, whereas the SMC has a more 
quiescent, smoother appearance. It seems surprising that the SMC dataset yields a so 
much larger shell catalog than the LMC. It turns out that for the LMC, SFR/SFRcrit ~ 1, 
whereas for the SMC, SFR/SFRcrit ~ 0.1 (Oey 2001). It wiU be interesting to examine 
shell catalogs and properties for larger samples of galaxies, from the THINGS H i survey 
(E. Brinks, these proceedings), for example. 

5. Detailed correspondence with star-forming regions 

Last, but not least, we seek a one-to-one correspondence between the observed super- 
bubbles and parent star clusters. To date, results are not as clear-cut as is desireable, 
but there does appear to be broad consistency, though somewhat controversial. Kim et 
al. (1999) find enough correspondence between the LMC Hi and Ha data to suggest an 
evolutionary sequence based on the morphological and kinematic relationship between 
these. Hatzidimitriou et al. (2005) examine the spatial relationship between the H i shells 
in the SMC and existing data for young stellar clusters, and do find broad correspon- 
dence, although the quantitative significance of the correlations is ambiguous. They also 
find that about 10% of the Hi shells show no counterparts in the stellar population. In 
Holmberg II, a targeted search for clusters in the H i shells failed to find the expected 
correspondence at optical wavelengths (Rhode et al. 1999), although Ha and FUV ob- 
servations seem more consistent with a mechanical feedback origin (Stewart et al. 2000). 

Finally, in the theme of this Symposium, triggered star-formation is also an association 
of massive star clusters with the creation of superbubbles. There are many examples of 
two-stage, sequential star formation; some well-known examples are Nil (Walborn & 
Parker 1992), N44 (Oey & Massey 1995), and N51 D (Oey & Smedley 1998; Cooper et al. 
2004) in the LMC, and the Rosette Nebula (Williams et al. 1995) in the Galaxy, among 
others. While it is difficult to establish a causal relationship between two-stage star- 
forming regions, a three-stage sequence is much more convincing. We recently identified 
the Perseus W3/4 complex in the Galaxy, identified as a triggered system by Thronson 
et al. (1985), as three-stage, hierarchical star formation, whose morphology is difficult to 
interpret as anything other than causal, triggered star formation (Oey et al. 2005). 

6. Summary 

To summarize, we see that observations are broadly consistent with mechanical feed- 
back from the most massive stars being responsible for the formation of most superbub- 
bles and shell structures seen in star-forming galaxies, and that these objects can be 
understood in terms of the standard, adiabatic evolution driven by massive star winds 
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and supernovae. There is empirical, multi- wavelength confirmation of the different tem- 
perature phases associated with the shock heating of superbubble interiors and also 
photoionization by OB stars in the youngest objects, as is qualitatively predicted by 
the adiabatic model. Statistical properties of shell populations are quantitatively consis- 
tent with predicted distributions in size and expansion velocity, based on known global 
properties of the parent star-forming regions and young star clusters. Above a thresh- 
old star-formation rate, we predict interaction among the shells and shredding of the 
ISM, which appears to be observed in the contrasting H l morphology and shell popula- 
tions between the LMC and SMC. The expected one-to-one correlation between young 
massive-star clusters and superbubbles remains somewhat ambiguous, although a general 
consistency is tentatively seen. 

When comparing the observed evolution and kinematics of individual objects with 
standard predictions based on detailed knowledge of the parent stellar population and 
other input parameters, we find that the objects are also broadly consistent with adi- 
abatic evolution, but that the shells are invariably smaller than expected. There are 
a number of possible reasons for this growth-rate discrepancy. (1) The assumed input 
power may be overestimated, especially with respect to measured mass-loss rates for 
stellar winds. (2) It is possible that the ambient density has been systematically under- 
estimated, although our resolved observations of H i environments for three objects does 
not especially support this interpretation. (3) There may be a systematic underestimate 
of the ambient pressure, since the various contributors to interstellar pressure all corre- 
late with star-formation activity. (4) Enhanced radiative cooling may be occurring in the 
hot superbubble interiors, caused by mass-loading or metallicity enhancements from the 
parent SNc. (5) Finally, a somewhat-overlooked mechanism is the transfer of mechanical 
energy from the superbubbles to cosmic rays, which is plausible since superbubbles are 
an especially effective acceleration environment. 
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